Scalar resonances in a unitary tttt S-wave model for t\^t[ tt"*" 

D. R. Boito/'^'* J.-P. Dedonder,^ B. El-Bennich,^'^' O. Leitner,^'^' * and B. Loiseau^ 

^ Grup de Fisica Tedrica and IFAE, Universitat Autdnoma de Barcelona, 
E-08193 Bellaterra (Barcelona), Spain. 
^Instituto de Fisica, Universidade de Sao Paulo, 
-Q^- CP. 66318, 05315-970, Sao Paulo, SP, Brazil. 

I ^Lahoratoire de Physique Nucleaire et de Hautes 

(N 

\ Groupe Theorie, Univ. P. & M. Curie, 

ly-^ ' ^ PI. Jussieu, F-75252 Paris, France 

"^Physics Division, Argonne National Laboratory, Argonne, IL, 60439, USA. 



(N 



i-S^ ' ^INFN, Laboratori Nazionali di Frascati, 

Via E. Fermi 40, 1-00044 Frascati, Italy 



(Dated: February 25, 2009) 



_ , Abstract 
CN ■ 

' We propose a model for D+ — > 7r+7r~7r+ decays following experimental results which indicate that the two- 

cn ; 

Tlj" , pion interaction in the 5- wave is dominated by the scalar resonances /o(600)/cr and /o(980). The weak decay 

00 ■ . 

■ amplitude for _D+ Rn^ , where i? is a resonance that subsequently decays into tt+tt , is constructed in a 

OA ' factorization approach. In the S'-wave, we implement the strong decay R — s- tt+tt^ by means of a scalar form 

QQ I factor. This provides a unitary description of the pion-pion interaction in the entire kinematically allowed 

mass range rn^^ from threshold to about 3 GeV^. In order to reproduce the experimental Dalitz plot for 
L)+ 7r+7r^7r+, we include contributions beyond the S'-wave. For the P-wave, dominated by the p(770)°, we 



o 

> 
X 

^ ' use a Brcit-Wigner description. Higher waves are accounted for by using the usual isobar prescription for the 



/2(1270) and p{lA50)'^. The major achievement is a good reproduction of the experimental m'^^ distribution, 
and of the partial as well as the total 7r+7r^7r+ branching ratios. Our values are generally smaller 

than the experimental ones. We discuss this shortcoming and, as a byproduct, we predict a value for the 
poorly known D a transition form factor at = ra^. 
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I. INTRODUCTION 



In 2001, the E791 collaboration found a very strong evidence for a light and broad scalar- 
isoscalar resonance in — 7r+7r~7r+ decays [1], confirming the existence of the /o(600) (also 
referred to as the a). This pioneering work was soon followed by the authentication of another 
elusive scalar, the i^Q(800) (or k), in D+ — K'tt'^it^ [2]. In the past, several analyses of tttt 
scattering data already claimed the presence of the a in the form of a pole close to threshold 
and with a large imaginary part [3-5]. However, its manifestation in scattering is subtle and 
one can consider the E791 experiment as the first solid empirical evidence for this resonance. The 
understanding of this pole in tttt scattering has improved considerably thanks to Chiral Perturbation 
Theory (ChPT) [6] and to dispersion relations, namely Roy's equations [7] (for recent works see, for 
instance, [8-10]). In 2006, the a pole was obtained from a theoretical analysis combining these two 
ingredients and yielding the accurate result = (4411*1^2 — z 272^^2 5) MeV [11]. Furthermore, 
in the last few years, experimental evidences for the a in other processes such as J/ijj ^ wtt+tt" 
[12] and, notably, from new analyses of vr^vr'vr^ decays [13, 14] have been published. Thus, 

at present, one can safely state that the a is the lightest resonance in the hadronic spectrum. 

However, in spite of all theoretical progress, a comprehensive description of the reaction — > 
7r+7r~7r+ is still to be accomplished. The reasons are twofold. Firstly, the c-quark mass lies in an 
intermediate range of energy, between the realm of light quarks (u, d and s) and that of heavy 
quarks {h and t). Although, on the one hand, decays of light mesons such as the kaon can be 
treated within the framework of ChPT and, on the other, decays of the B can be calculated within 
QCD Factorization [15], Heavy Quark Effective Theory [16] and Soft Collinear Effective Theory 
[17], no such rigorous framework exists for the D. Secondly, one deals with a three-body final state 
which renders a full treatment of mesonic final state interactions (FSI) most involved. 

In the search of a more sound theoretical framework to treat this type of reaction, experimen- 
talists usually fit the Dalitz plot for the three-body decay with the isobar model. Schematically, 
the model consists of a trial amplitude of the form 



where the first term in the r.h.s. corresponds to a non-resonant background and the sum runs over 
all the n resonances that contribute to the decay. In Eq. (1), the parameters and (f)i are real 
constants and the sub-amplitudes Ai (depending on invariant masses) whose main ingredients are 
relativistic Breit-Wigner functions (BW), represent the propagation and decay of each resonance. 
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In the course of the analysis, E791 discovered that the usual set of tttt resonances was not sufficient 
to yield a good fit and, therefore, the addition of a new BW was required. Leaving its mass 
and width free to float, an improvement of the least square function, x^, was obtained with the 
parameter values = (478+23 ^ 1^) MeV and = (324^4^ ± 21) MeV. These values, however, 
depend strongly on the specific BW chosen to fit the Dalitz plot and, in the literature, one finds 
various choices of BW propagators. Moreover, this simple description for a very broad resonance 
close to a threshold, like the a, has many deficiencies that are discussed, for instance, in Ref. [18]. 

On the theoretical side, a possible way of improving Eq. (1) consists in replacing the BW func- 
tions with expressions based on the knowledge of scattering amplitudes. This procedure respects 
unitarity and reveals the relation between scattering and production experiments. Using a descrip- 
tion of the /S-wave two-body FSI within a unitarized ChPT framework [19], Oiler has proposed, 
in Ref. [18], a modified version of Eq. (1). Albeit successful, this model does not tackle the weak 
vertex of the reaction and, consequently, the constants and (pi in Eq. (1) remain fit parameters. 

Motivated by the study of CP violation, more detailed versions of Eq. (1) were produced in 
the context of three-body hadronic B decays, in which both the weak and strong interactions are 
treated: B tttttt was considered in Refs. [20, 21] whereas B K-rtt and B — KKK were treated 
in Refs. [22-24]. In these works, the weak decay is evaluated with the help of the effective weak 
hamiltonian within QCD factorization whereas the hadronic FSI are taken into account by means of 
unitary Kir and tttt form factors constrained by scattering data and ChPT. These models confront 
data very well and are precisely the basis of the approach we follow in the present work. Our main 
purpose is to test the description of tt+tt" pairs in a relative 5- wave state in vr+vr^vr"'". 

The decay amplitude for — >■ 7r"'"7r~7r"'" is constructed as follows. We assume that the three- 
body decay is always mediated by a resonance R as suggested by experimental analyses [1, 13, 14]. 
Then, for the tt+tt" pairs in an S-wave state, denoted hereafter (7r+7r'")5, we factorize the decay 
amplitude Rtt'^ using the effective weak hamiltonian within naive factorization. Afterwards, 

the three-body (7r+7r~)5 7r+ final state is constructed from the intermediate i?7r+ state employing 
the TTTT scalar form factors introduced in Ref. [22]. This ensures that our description of FSI is 
unitary and includes the coupling to the KK state. The form factor is based on the experimental 
scattering phase shifts (5,r7r and (5^^ previously studied in Ref. [5]. As usual, we work within the 
quasi two-body approximation in which interactions of the remaining tt"*" with the (7r"'"7r~)5 pair 
are neglected. This procedure is repeated for tt+tt^ in a P-wave, which is well approximated by 
the /9(770)° resonance. Finally, two higher mass resonances, the /2(1270) and the /9(1450)°, are 
included phenomenologically using the isobar model. Our final amplitude is then fitted to the 
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signal function employed by the E791 collaboration [1]; this comparison is carried out following a 
scheme presented in Ref. [18]. 

The present paper completes the description briefly reported in Ref. [25] and is organized as 
follows. In section II, we discuss the effective weak hamiltonian and the weak amplitudes employed 
in our description for the S- and P-waves. The construction of the three-body final state is pre- 
sented in section III, results for our fits are displayed in section IV and a summary and discussions 
are given in section V. 

II. WEAK AMPLITUDES 

Our phenomenological description of weak decays involving S- or P-wave resonances is based 
on the effective weak hamiltonian, H^g, which is obtained by integrating out the heavy degrees 
of freedom of the Standard Model (SM) lagrangian. The hamiltonian is written as an Operator 
Product Expansion (OPE) and reads 



where Gp = 1.16637(1) x 10 ^ GeV ^ [26] is the Fermi decay constant, Vckm are products of 



Furthermore, /x is the renormalization scale which, in our case, is taken to be of order rric- The 
operators Oi{^) represent the local four-quark weak interaction in the effective theory while Ci{ji) 
describe the hard or short-distance physics and arc calculated pcrturbatively in the full theory, in 
this case the SM. The contributions of all particles with mass m > /i = nic, such as the heavier b 
and t quarks and the W bosons, are included in Ci{fi). 

In the present case, we need to consider the Cabibbo suppressed transition c dud. Conse- 
quently, the tree level matrix elements are governed by the coefficient V^^^ = VcdVud = ^d, which, 
in Wolfenstein's parametrization [27], is of 0{X) where A = 0.2257 [26]. In principle, our amplitudes 
receive contributions from strong penguins as well but, using the unitarity of the CKM matrix, 
one sees that these are governed by V^^.j = VcbV*f^ which is of 0{X^). Thus, penguin amplitudes 
are strongly CKM suppressed and can safely be neglected. This suppression contrasts with the 
situation found in analogous B decays where, since one explores a different sector of the CKM 
matrix, penguin operators give rise to sizable contributions [20]. 

Taking into account only tree operators, the matrix element of H^ff for the decay i?7r+ 




(2) 



CKM matrix elements, Ci{^) are Wilson coefficients and Oi(/x) local operators entering the OPE. 
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can be written as 

2 

{R'K+\n,Q\D+) = ^\dY.'^,{^i){R'K+\^i\D+){^Ji), (3) 

where the operators 0\ and O2 are given by 

01 = ^71.(1 - 75)c U7''(l - 75)d, 

02 = u-l.{l - 75)c (i7"(l - l^)d. (4) 

To deal with the matrix element (i27r^|Oj of Eq. (3), we assume that factorization at 

leading order (in A.QCD/fnc and a^) holds 

{R'K+\n^^\D+) = ^Arf[ai(/i)(7r+|iZ7,(l -75)d|0)(i?|dV(l -75)c|I)+) 

+ a2{n){R\d-f,{l - 75)d|0)(7r+|u7'^(l - 75)c|£'+)] . (5) 

In this context, the new coefficients ai{ji) and 02 (/u) which arise are expressed in terms of the 
Wilson coefficients Ci {/j,) and C2 (/x) and of the number of colors iVc = 3 as 

ai(/x) = Ci^i) + ^C'2(//), a2(/x) = C2(m) + ]^C'i(^). (6) 

We then have products of non perturbative hadronic matrix elements that assume a decomposition 
in terms of Lorentz invariant form factors. At this point, a remark is in order. Since rric is smaller 
than nif, by roughly a factor three, non-factorizable contributions of order Kqcd/tt^c are larger 
than in B decays. Hence, the factorization approximation may be less reliable for D physics. 
Factorization, however, was applied successfully to two-body D decays in the seminal papers by 
Bauer, Stech and Wirbel [28]. Later, it was used in the 80's [29] and 90's [30] to describe three-body 
hadronic D decays and, more recently, in 2001, Dib and Rosenfeld [31] have used the factorization 
approximation to treat the air'^ decay in order to obtain, from the then novel E791 data, 

the CTTTTT coupling and the form factor for the transition D ^ a. Finally, the transition D 
/o(980) was recently studied employing the same factorization scheme [32]. Thus, in spite of the 
complications introduced by the c mass scale, this approach is a reasonable starting point for a 
phenomcnological analysis. Moreover, it enables the description to benefit from the treatment of 
mesonic FSI developed in the context of B decays^ . 



^ This is the reason why we do not adopt the effective mesonic lagrangian of Refs. [33, 34]. 
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A. D+ cnT+ and D+ /o(980)7r+ decays 



Prom the available experimental analyses it has become clear that the tttt S'-wave gives the 
most important contribution to the decay 7r+7r~"7r+. The a was found to account for almost 

50% of the decays whereas the /o(980) accounts only for about 6% [1]. This picture remains 
unchanged with the more recent analyses of FOCUS [13] and CLEO [14] (see Table III for the 
values extracted from these papers). In addition, there are indications for a sizable component 
arising from the higher mass scalar resonances /o(1370) and /o(1500), although the data analyses 
are not conclusive. The E791 analysis includes only one state, the /o(1370), whereas the CLEO 
collaboration claims that both are necessary to yield a good fit within the isobar model. In our 
model, the a and the /o(980) enter explicitly in the amplitude while one higher mass state close to 
1500 MeV is accounted for as a pole in the (7r"'"7r~)5 form factor, discussed in section III. 

The exact nature of scalar mesons, i.e. whether their wave function is dominated by qq or 
(fq'^ states, or even glueballs, is not yet elucidated. In the case of the a, the situation is even 
more obscure, since the issue whether it is a pre-existing quark state or dynamically generated 
by the strong S-wave tttt interactions is still under debate. Therefore, the inclusion of the a 
in our amplitude in far from obvious. On the other hand, we know that the a is an isoscalar 
state and no indications of a strange quark component in its wave function exists. Thus, we 
assign to the a the minimal quark content compatible with these two facts, namely {uu + dd)/\/2. 
Concerning the quark content of the /o(980), it is well established that it has a non-negligeable 
strange component. Higher Fock states, e.g. q^(^ and cf'q^g, may also be important to achieve a 
comprehensive description of its wave function^. In the context of heavy meson decays, however, 
there are indications that the qq component should be dominant [35]. Here, we consider the /o(980) 
as pure qq while allowing for an admixture of ss: 

fo ^ + ^^"f ^™ ^mix + ss cos ^mix, (7) 

where ^niix is the mixing angle. With these definitions, we can compute the weak amplitude for 
D+ cr7r+ and D+ /o(980)7r+ decays. 

The matrix elements (i?|g7ju(l — 75)^10) where R = a ov /o(980) and q = u,d ot s vanish by 
C invariance. Hence, the decay amplitude has no a2 contribution (see Eq. (5)). The annihilation 
topology is likely to be neglected since it contains form factors of light mesons evaluated at high 
momentum {q^ = mj^) [31]. The weak amplitude is then purely proportional to ai. Finally, we 

^ For a more complete discussion see Ref. [32] . 
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note that, in the next section, we construct the three-body final state from the intermediate air'^ 
and foTT'^ states with the help of the tttt scalar form factor, thereby properly taking into account 
the strong FSI in this channel. 

The relevant form factor for the transition D ^ a, denoted by F^^'^{q'^), is defined as^, 

q^{a{p)\d^,{l - ^,)c\D+{Pd)) = -i{ml - ml)F^-^%q% (8) 

where = {PD—pf and F^^^^q"^) = FQ^'^{q^)/\/2. Furthermore, with the usual decomposition 
one has 

{'K+{p)\uYi}-lM^)=iUl^- (9) 

It is then straightforward, once factorization has been applied, to write the amplitude for — 
(TTT"'" with the help of Eq. (5) 

A{D ^ a7r+) = (a7r+|Weffl^+> = ^A.aiM {ml - ml) F^--{ml). (10) 

This expression coincides with the one found in Ref. [31] and can easily be compared with the 
amplitude for B ^ an from Ref. [20] as well. However, one should note that in Refs. [20, 31], the 
normalization factor 1 / \/2 is not included in the a wave function. Finally, we have computed the 
strong penguin contributions to Eq. (10) and we have checked explicitly that they can be neglected 
due to CKM suppression. 

For the decay — /o (980) 7r+ we obtain the following result 



A{D+ - /o7r+) = ^ Arf ai(me) U {ml - m%) Pl^^^^ml), (11) 
where we have defined 



Fr''^{ml) = '^Fr^\ml). (12) 
B. p(770)° 7r+ decay 

The evaluation of the amplitude A{D'^ — /9(770)°7r+) from Eq. (5) is quite analogous to the 
previous ones except that the contribution of the colour suppressed tree diagram does not vanish. 
The quark content in this case is = {uu — dd)/V2 (for simplicity we denote by the ^(770)^^). 



^ Throughout this paper, we employ the decompositions and definitions of Refs. [20, 21] for the hadronic form 
factors, unless otherwise stated. 
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In the colour allowed term we need to consider the transition D p which can be parametrized in 
terms of form factors using the general P amplitude, where P and V represent a pseudo-scalar 
and a vector meson respectively: 



(13) 



Here q = pp — pv, is the vector meson polarization. Furthermore, and represent bilinear 
vector and axial-vector quark currents, respectively. In the colour suppressed topology, we need 
the transition D ^ tt obtained from the general Pi P2 transition 



{P2{p'pM - j^lPiipp)) 
We define the vector decay constant as 



, , , — rnk 
{Pp+Pp)h 2 



2 2 
rnr — m 



{V\j^\0)=mvfve;, 



F^iq') + ^^^^q,F,iq'). (14) 



(15) 



where we have included the factor my to have a decay constant fv with dimension of energy^. 
With these definitions the result reads 



^(L>+^ A+)= -2(e*-g)r?o, 



(16) 



where 



\ 
1 = -TT^dmp 



(17) 



and q = {Pd —Pp) = P-n- The minus sign in Eq. (16) arises from the fact that the dd component of 
the is the one that intervenes. This result has the same structure as the one found in Ref. [20] 
for the B~ p^iT~ decay. 



III. HADRONIC FINAL STATE INTERACTIONS 



In this section wc describe the construction of the thrcc-pion final state from the weak amplitudes 
Eqs. (10), (11) and (16). This is done in order to take into account the FSI and the three-body 
phase space. By FSI, we mean the mesonic interactions in the final state after hadronization. 

This last definition is not the same as in Refs. [20, 21]. 
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It is important to define the kinematics we employ. We are considering the generic amphtude 
M = M{D'^ — Tr+TT-Tr"*") with four-momenta labeled as follows 

D+{pd) ^ 7r+(pi)7r-(p_)7r+(p2). (18) 

Since in the final state we have two identical vr^, the amplitude has to be symmetric under 
the exchange pi ^ P2- What is more, in order to work with a Lorentz invariant Daliz plot, it is 
convenient to define three invariant combinations of momenta 

s = (^^l+P2)^ t = {p-+pif, u = {p-+p2f. (19) 

They correspond to the usual Mandelstam variables with 

s + t + u = ml, + ^ml±, (20) 

so that only two of them are independent. Resonances occur only in the t- and -u-channcls since 
there are no isospin 2 resonances. We denote by Ar{u, t) the amplitude for a decay mediated by a 
resonance R in the u-channel. The final symmetric result M.R{u,t) is hence obtained by summing 

MR{u,t)=AR{u,t) + AR{t,u). (21) 

A. S-wave 

The FSI are taken into account by means of the tttt scalar form factor. This method was 
introduced in Ref. [36] and was later applied to B ^ tttttt decays in the vicinity of the a pole 
in Ref. [20] using the form factors obtained in the context of unitarized ChPT [19]. In Refs. 
[22, 23] a similar description was used to describe the /o(980) m. B ^ ttttK and B KKK 
employing a different set of form factors that rely on a previous analysis of tttt and KK scattering 
data [5]. These form factors are unitary and contain both the tttt and KK channels. Here, we 
briefly summarize the model since all the details can be found in Ref. [20] and in the appendix of 
Ref. [23]. 

The two pion scalar form factor, r"^(a;), that is relevant to our work is defined as [36] 

(0|nn|7r(p)7r(y)) = V2Bor"(x), (22) 

where nn = (uu + dd)/\/2, Bq is proportional to quark condensate Bq = —{0\qq\0)/f^ and x = 
{p+p'Y- Since we want to describe — 7r+7r~7r+ and we have defined in Eq. (10) the amplitude 
for the transition D'^ — ^ u-k'^ , we need to introduce a function, Jl^jT^^ix), that describes the a 
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propagation and decay^, i.e. the final state interactions. The full amplitude air'^ 
(7r+7r~)5 7r+ is given by 

Aa{u,t) = A{D+ aTr+)Uan^{u). 

The a {it^tt^)s decay is described without resorting to BW expressions. It can be obtained 
from the complex conjugate of Eq. (22) assuming that, close to the a pole, this resonance gives the 
dominant contribution to r"(x). We have [20] 

The normalization constant 

V 3 {a\nn\0) 

is, in principle, unknown as it depends on the matrix clement {a\nn\0). 

To further clarify the meaning of IlaTr-K{x), it is convenient to consider its analogue in a BW 
framework 

where gairn is the coupling between the a and the pions. We can obtain an estimate for Xa 
comparing expressions (23) and (25) at a; = m^. We obtain 

We take the central values rrifj = 478 MeV and To- = 324 MeV from the E791 fit and use gann = 
2.52 GeV [37]. This yields, with the form factor of Ref. [22], 

Xa ~ 29 GeV~^. (27) 

Prom the weak amplitude given in Eq. (10) and from the expression of no.7r7r(a;) the amplitude 
for Z)+ — (T7r+ — (7r+7r~)5 7r'*" reads 

AAu,t) = ^Adai(me)/.(m^+-m2)Fo^--(m2)x.r"*(u), (28) 

where the t dependence is implicit. Similarly for the intermediate resonance /o(980), one has 



Af,{u,t) = ^\da,{rn,)U{rnl+-ml)F^-^f''{ml)xfor^*{u). (29) 



nCT7r7r(a;) corresponds to TaTrTr{x) of Ref . [20]. 
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For the weak decay amplitude we make use of Eq. (11) and one should note that the normalization, 
denoted Xfo^ which by virtue of Eq. (26) is proportional to the coupling 5/g7r7r, differs from the one 
found in Eq. (28). 

Thus far, we have considered only a small energy range around the resonance poles. In other 
works, this prescription for the 5'-wave was always used only within a limited region of the spec- 
trum [20, 22]. We want, however, to describe the whole Dalitz plot for D'^ tt~^tt~tt~^ where the 
TTTT invariant mass range from 2m^± < ^/s < {rnjj+ — m^±), i.e. between 280 — 1700 MeV. To 
this aim, an ansatz is required to provide us with an expression for the entire S-wave. Since the 
amplitudes (28) and (29) are both proportional to r"*(n), we propose the following amplitude for 
the S'-wave 

As{u,t) = --j=Xdai{mc) U {^0+ - u) Xeffr"*(u). (30) 

In the last equation, XcK is a new normalization constant that encompasses all the form factors and 
normalizations for the scalar resonances. In addition, we have replaced the terms (m^+ — m|j) by 
(m|,+ — u) ; this -u-dependence suppresses the contribution of higher mass resonances [22] . It is not 
easy to obtain a good estimate for Xefi since it receives contributions from all the scalar-isoscalar 
states. The following lower bound of XeS results from Eqs. (28), (29) and (30) 

sin ^jnix TjiD^/o / 2 \ 
Xeff > ^ -^0 (^tt) X/o ■ 

Using F(f^-^°(m2) = 0.215, average value from the two models of Ref. [32], and Xfo = 28.9 
GeV~^ [23] one then gets the estimate 

Xeff > 6.2 GeV-^ (31) 

Prom Eqs. (30) and (21) we can construct our final expression for the — (Tr+Tr")^ 7r+ 
amplitude 

Msiu, t) = ai(mc) U Xeff [{ml+ ' u)T^*{u) + (m^+ - t)V^\t)\ . (32) 

The last expression has only one parameter that can be considered as unknown: Xefi- ^ 
r"'*(x) is concerned, with x = u,t, we employ the form factor of Ref. [22]. It is unitary and takes 
into account the coupling to the KK channel^. This form factor is obtained within an on-shell 



In a two coupled channel (tttt, KK) description of the final state interactions for decays, the FSI are incorpo- 
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approximation and can be written explicitly in terms of the S-wave scattering phase shifts (57r7r(a;) 
and Sf(^{x) as 



(33) 



where 'q(x) is the inelasticity for tttt scattering, ki = s/x/A — and k2 = — m|.. For 5.j^t^{x), 

Sxx{x) and rf{x), we employ the results of Ref. [5]. In addition, Eq. (33) depends on the production 
functions i?5J^(x) and R^j^{x) introduced in Ref. [36], which result from a matching to the ChPT 
expansion of r"(a;). In practice, these functions can simply be written as (a;) = Cj + diX. The 
real coefficients q and dj, which depend on the low energy constants of ChPT, were determined in 
Ref. [36] and updated in Ref. [38]. However, the validity of the production functions beyond ~ 1.2 
GeV is not guaranteed. To circumvent this problem, we introduce in our form factor a cut-off Xcut 
above which we saturate the Rf{x), namely 

{Ci + diX for X < Xcut, 
(34) 
Ci + diXcut for X > Xcut- 

This procedure does not affect the unitarity of the form factor but, of course, it introduces an 
additional parameter Xcut in the model. Our fits are done for different values of Xcut in order to 
carefully ascertain the dependence on this parameter. 



B. P-wave 



The construction of the three-pion final state from the intermediate p^Tr'^ state is similar to that 
done for the S-wsve. However, in the case of the it is not crucial to employ the vector form factor 
of the pion to describe the FSI. Since the p^ is a relatively narrow resonance, far from threshold, that 
strongly dominates the corresponding form factor, the BW description gives a good approximation. 
Prom the coupling of the /9° to the pair tt+tt" defined as (7r+(g+)7r~(g_)|p°) = gp{q- — q+) [20] 



rated in the form factor via the following unitary equation system (Eq. (11) of Ref. [22]): 

rr{x) =R?{x)+ j2{ki\R]{x)Gi{x)Tii{x)%), 

where \ki) and \kj) represent the wave functions of two mesons in the momentum space and the indices i,j = 1,2 
refer to the tttt and KK channels, respectively. The matrix T is the two-body scattering matrix and the functions 

Gj{x) arc the free Green's functions. With this definition, the form factor of Eq. (33) corresponds to ri*{x). The 
driving terms entering in these equations are given by production functions R"{x) representing the meson-meson 
formation from qq pairs. Further details can be found in Ref. [22]. 
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and with the use of Eq. (16) we have 

Apo{u, t) = r]o{t- s) np^7r(w), (35) 

where rjo is given in Eq. (17) and s, t and u are defined in Eq. (19). The factor {t — s) comes from 
the sum over the polarizations of the p° and the function HpTTTr {u) is defined by 

np..(«) = — (36) 

nip — u — impTp(u) 
For the running width Vp{u) we take the usual relativistic prescription 



W-^rr^^--jJ , (37, 

where p{u) = \/u/A — and r^°* is the total decay width. The final expression for the amplitude 
of the decay mediated by the is then 

Mpo {u, t) = rio [{t - s) n^^(n) + {u-s) n^^(t)] . (38) 



IV. RESULTS 

A. Parameter values 

Since the experimental situation of the P-wave in — p° tt"^ is less controversial than that 
for the S-wave, this channel can be used to ascertain the quality of the model. With Eq. (38) we 
can calculate the branching ratio for the p° 7r+ decay and compare it to the experimental 

average [26] 

^PDGp+ ^ pO^+^^o ^ ^+^-^ ^ (g_2 ± 1.5) X 10-1 

For the numerical input needed we take: ai{mc) = 1.15 [39], a2{mc) = —0.25 [39], = 
130.4 MeV [26], fp = 0.209 GeV [23], ^/^^^(m^) = 0.8 [40], A^-^^iml) « ^^^''(O) = 0.75 [41], 
Qp = 5.8 [20]. Using the standard formula for the three-body decay rate [26] and taking into 
account the symmetry factor 1/2 we have, using Eq. (38), 

H(L>+ ^ p°7r+, p° ^ TT+TT-) = 8.63 X 10"! (39) 

Our branching ratio is rather sensitive to the value of the form factor AQ^^im^) which has an 
uncertainty of about 20% at = [41] . Therefore, the theoretical error associated with our result 
is large. However, since the central values for the parameters yield a result in agreement with the 
experimental average, we keep these values. Thus, the amplitude M.po{u,t) serves as a benchmark 
to the determination of the other parameters of our model. 
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B. Fits to E791 signal function 



Since we do not have the real data at our disposal, the fit procedure consists in reproducing the 
E791 signal function and comparing our model to it. To this aim, we closely follow the method 
of Ref. [18]. Aitala et. al. in Ref. [1] used in their best fit the trial amplitude Eq. (1) with six 
resonances, namely the a, the /o(980) and /o(1370) whose quantum numbers are l'-^{J^') = 0+(0"'"), 
the p(770)° and p(1450)° with 1+(1") and, finally the /2(1270) with 0+(2+). Let us remind that 
the amplitude has also a complex constant ajv^e*'^^^ identified with the non-resonant background. 
With these ingredients, the signal function can be written 

6 

M^'^'iu, t) = aNRe'^^^ + ^ aie'^^ [Af''{u, t) + Af'^it, u)] , (40) 

i=l 

where the individual amplitudes Af^^^iUjt) are modeled as 

Ar\u,t) = F{){u) X Fi{u) X ^i{u,t) X BWi{u). (41) 

In the last expression, F^{u) and F/{u) are Blatt-Wcisskopf damping factors that depend on the 
spin J of the resonance, il^ (n, t) are angular factors and BWi{u) are the Brcit-Wigner propagators. 
The full expressions for these functions are rather lengthy and can be found in the original paper [1] 
or, more detailed, in Ref. [18]. Prom Eq. (40), and employing the fit results for aj and [1] as well 
as the values for masses and widths used by the E791 collaboration we are able to reproduce the 
signal function. Then, we generate a Dalitz plot with points separated by 0.05 GeV and normalize 
this plot to the number of observed signal events, 1124. Since the Dalitz plot is symmetric under 
u ^ t it is sufficient to fit only half of the plot. 

Our complete amplitude for the S'-wave is given by Eq. (32) and the decay mediated by the 
is described by Eq. (38). With only these two contributions, however, it is not possible to achieve a 
reasonable reproduction of the E791 signal function. We have to include the two other resonances 
that give sizable contributions to the fit: the /2(1270) and the p(1450) (denoted for simplicity /2 
and p' respectively). This is done with the help of the isobar model, using the same expressions 
as those of [1]. We leave free, in our fit, the corresponding magnitudes a/j and api and phases 
and (f)pi. The final amplitude for the decay -D+ — 7r+7r~7r+ reads then 

M{u,t) = Msiu,t) + Mpo{u,t) + Mp'iu,t) + Mf^{u,t), (42) 

where Msiu,t) is given by Eq. (32) and Mpo{u,t) by Eq. (38). For Mp'{u,t) we use 

Mp'{u,t) = ap^e't'p' [A^/'{u,t) + A^pJ''{t,u)] , (43) 
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where the explicit form of Ap' is given in Eq. (41). Analogously, for Mf^iujt) we use 

Mf,iu,t) = af^e'-^h [Ajl'\u,t)+A%'\t,u)] . (44) 

We shall refer to this model as "model A" . 

The amplitude Eq. (42) receives contributions from five different resonances. These arc: the 
three scalars a, /o(980) and /o(1500) that appear as poles in our 5-wave form factor Eq. (33) [22], 
the P-wave resonances p(770) and p(1450) and the D-wave represented by the /2(1270). Note 
that we do not include a non-resonant (NR) amplitude. The necessity for such an amplitude is 
controversial. The E791 collaboration has shown that the inclusion of the a reduces the contribution 
of the NR background to less than 10%. More recently, the CLEO collaboration did not find any 
significant evidence for the NR amplitude and an upper limit of 3.5% was established [14]. In 
addition, the NR amplitude can be energy dependent and the simple complex constant evenly 
spread over the whole phase space may be an Tinreliable model. Most importantly, the vrvr scalar 
form factor T"* (x) already includes the NR contributions to nir scattering which could generate a 
double counting of the background. 

TABLE I: Fits of model A, Eq. (42), to E791 signal function Eq. (40). Uncertainties are solely statistical. 



y/¥^t = 1.0 GeV ^/i^ = 1.2 GeV = 1-4 GeV 



Xeff [GeV-i] 


6.5 ±0.3 


6.1 ±0.3 


6.0± 0.3 


X 10^ 


(4.7 ±0.5) 


(4.2 ± 0.4) 


(4.2 ± 0.4) 


0/2 (rd) 


-6.03 ±0.18 


-6.11 ±0.20 


-6.18 ±0.20 


a'p X 10^ 


(2.2 ±0.8) 


(2.6 ±0.7) 


(3.0 ±0.6) 


0P' (rd) 


-0.57 ±0.27 


-0.30 ±0.20 


-0.31 ±0.17 


xVd.o.f. 


0.20 


0.22 


0.22 



Only five free parameters occur in Eq. (42). Four of them (a/j, 4'f2i '^p' 0p') arise from the 
isobar model description of the /2(1270) and p(1450). The fifth is the real constant XcSi introduced 
in Eq. (30). Therefore, the relative weak phase of 5'- wave with respect to the (P is fixed, as well 
as the strong phases. In Table I wc display the results for fits to the E791 signal function. For the 
functions E!^^{x) and BJ^^{x) of Eq. (33) we use the updated values obtained in Ref. [38]. We 
show the results for three fits in which we vary the value of the cut-off Xcut introduced in Eq. (34). 
The fit parameters do not depend much on x^ut- When quoting final values, we take the fit with 
■s/Xcut = 1.2 GeV and include an uncertainty due to the dependence on this cut-off. In Fig. 1 we 
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show both the E791 signal function and the result of the fit with ^Xcut = 1-2 GeV displayed as 
Dalitz plots. The projection of these two functions is compared in Fig. 2. 
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FIG. 1: Dalitz plot representation of the signal function employed by E791, Eq. (40), (l.h.s), same repre- 
sentation for the fitted amplitude given in Eq. (42) with the parameters of Table I for ^Xcut = 1-2 GeV, 
(r.h.s). 
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FIG. 2: Projection of the signal function of the E791 collaboration (full dots) and the result of the fit (model 
A) for the parameters given in Table I with ^x^ut = 1-2 GeV. 



Concerning the parameters of the fit, from Table I we see that the normalization constant 
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for the S'-wave, XeS^ is well determined and confirms our expectations derived in Eq. (31). The 
magnitudes of the /2(1270) and of the p(1450) are well constrained by the fit as well. On the other 
hand, the phases for these higher mass resonances are not well determined. Finally, since M. {u, t) 
and Mpi{u,t) in our amplitude Eq. (42) are exactly the same as in the function we are fitting to, 
namely Eq. (40), the interpretation of the x^/d.o.f as a measurement of the quality of the fit is not 
reliable. 

From Figs. 1 and 2, one sees that the main discrepancy in the fit comes from the a region. 
This is most probably due to the off-shell effects that are not included in the form factor given 
by Eq. (33). The omission of these effects in no-7r7r(a;) may lead to an underestimation of the a 
peak [42]. The model gives a much better description of the /9° and the /o(980) peaks, both very 
prominent in the projection. In the higher energy region, where the /2(1270) the p{lA50) and 
higher mass scalar states are present, our description is reasonable, although not as good as in the 

and /o(980) region. 

C. D — > cr treinsition form factor 

It is desirable to extract the form factors for the transition I?+ —>■ cr7r+ and /o(980)7r+ 
from the fit. The continuous description of the S-wave, however, renders this extraction difficult 
since the form factors are embedded in XeS- Nevertheless, we advance here a model which aims at 
determining F^~*°'(m^) and F^'^-^°(rn^). We then need to separate the a and /o(980) contributions 
in the S-wave. The drawback is however a non-unified description of (7r"'"7r~)5. 

For energies within the elastic domain, x < 4m|', the form factor of Eq. (33) is proportional 
to cos57r7r(x) [22]. Thus, it has a zero at the point xq where St^t^^xq) = tt/2. Numerically, from 
the analysis of Ref. [5], we have xq ~ (0.828GeV)^. We split our S-wave amplitude at this point 
and consider that for energies x < xq the {Tr~^n~)s is dominated by the transition D a and 
hence described by Eq. (28), whereas for x > xq the dominant transition is D ^ /o(980) given 
by Eq. (29). In practice, we substitute Eq. (30) by Eq. (28) when x < xo and by Eq. (29) when 
X > Xq. We shall refer to this modified version of the model as "model B". 

With this modification, we can now consider the products Xcr-^(P~"^('^^) and Xfo^Q^^° i^V) 
from Eqs. (28) and (29) as free parameters of the fit. The results of the fit are shown in Table 
II. The calculated uncertainties take into account the three possible sources: statistics, chang- 
ing the set of Rf{x) and varying Xcut- The dominant error in the case of Xo-F(f ^'^(m^) and 
XfoF^^^^irn^) comes from the two different sets of R^{x) functions that we have at our disposal 
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from Refs. [36] and [38]. This procedure is conservative, since it probably yields an overestimation 
of the theoretical error. The central values are obtained from the most recent determination of 
Ref. [38]. In this fit, statistical uncertainties are larger due to the additional parameter of the 
model indicating that, within the present experimental constraints, the model B is somehow over 
parametrized as compared to model A. Since the plots for this model are similar to Figs. 1 and 2 
we refrain from displaying them here. 

TABLE II: Fit of model B with ^Jx^t = 1-2 GeV to E791 signal function Eq. (40). Uncertainties include 
all the possible sources (see text). 





(8.4 ±1.4) GeV-i 




(5.6 ±1.9) GeV-i 


af^x 10^ 


5.3 ± 0.8 


(rd) 


-6.5 ±0.4 


ap> X 10*5 


2.1 ±0.9 


0p/ (rd) 


-0.2 ±0.6 


\2/(l.().f. 


0.21 



Let us now use the results of Table II to obtain the transition form factor for D ^ a. In order 
to disentangle the product X(jFq~*'^ {t^V) calculate Xa from Eq. (26). Gardner and MeiBner 
gave Xa- = 20 GeV"-*^ in Ref. [20] while employing the production functions R^{x) from Ref. [38], 
we get X(7 ^ 22 GeV~^. Prom Eq. (26) one sees that Xa depends on quantities that are not well 
known such as Fa- and ma and we shall take these values as an indication. Nevertheless, with our 
value one obtains 

Fq^" {ml) = 0.38 ± 0.06. (45) 

This result is to be compared to 0.79±0.15 from Ref. [31], 0.57±0.09 from Ref. [43] and 0.42±0.05 
from Ref. [44] which is however evaluated at = 0. Our value is compatible with the last one, 
but smaller than those of Refs. [31, 43]. More detailed discussions on the form factor ^(f^'^(m^) 
are beyond the scope of this work. 

Concerning the /o(980), we use the recent determination of FQ^^^\m1) from the phenomeno- 
logical analysis of Ref. [32] in order to estimate Xfo- The latter can be compared with the val- 
ues from Refs. [22, 23] to check the consistency of both approaches. Prom Ref. [32] we have 
F^~^-^°(m^) = 0.215 and with the result for the product XhFQ^^°{fn1) obtained in the fit, this 
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yields 

X/o = 26±9 GeV-\ (46) 

This result agrees within uncertainties witli tlic ones employed in B decays: x/o = 23.5 GeV~^ 
and x/o = 33.5 GeV'^ in the two models of Ref. [22] and Xfo = 28.9 GeV-^ in Ref. [23]. This 
comparison is to be done with care since one always has the product of the normalization x/o and 
the transition form factors. Therefore, smaller values for x/o can be compensated by larger form 
factors. 



V. SUMMARY AND DISCUSSIONS 



It is elucidative to quantify the weight of the different contributions in our model and compare 
them with other results in the literature. This is usually done through the fit fractions, defined, 
for a given contribution R as, 

f^dudt\MR{u,t)\'^ 
Ududt\EiMi{u,t)\''' ^ ' 

where T> indicates that the integrals are to be performed over the whole Dalitz plot. The fit 
fractions give an idea of how important a given resonance or partial wave is to the total decay 
amplitude. The sum of the fit fractions does not necessarily add up to one due to interference 
effects. However, one expects that this sum should not deviate widely from unity. In Table III we 
compare the fit fractions of models A and B with those of other works as well as with the values 
quoted by the PDG [26]. In our fit fractions we included the uncertainty from y/Xcut as well as 
that arising from the use of different sets of Rf{x) from Refs. [36] and [38] 

To clarify the content of Table III, one should remark that the E791 collaboration employed 
the isobar model [1], as already commented. Note that from the CLEO analyses [14], we only 
quote the isobar model results. FOCUS has performed an analysis where the S'-wave is described 
continuously by means of a i^-matrix previously determined from vrvr scattering data [13]. This 
model, therefore, is the closest to ours since both are based on tttt scattering analyses and share the 
on-shell approximation. Oiler's results [18], based on a description of FSI obtained in the context 
of unitarized ChPT [19], also rely on previous analysis of tttt scattering and include off-shell effects. 
In his model, the non resonant amplitude is kept and, although the resemblance to E791 data is 
striking, the price to pay is a huge interference that makes the sum of fit fractions to be 186% 
with 108% of TTTT S'-wave in the final state. Finally, we do not sum the fit fractions from the PDG 
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TABLE III: Fit fractions (in %) from E791 [1] , FOCUS [13], CLEO (isobar model) [14], Oiler [18], PDG [26] 
and models A and B. The uncertainties have been summed quadratically. Results marked with an asterisk 
are the sum of all S-wave contributions and are extrapolated from the original works. From Ref. [18] we 
included the 6% of KK — > (7r7r)s in the (TrTr)^ fit fraction. 





E791 


FOCUS 


CLEO 


Oiler 


PDG 


mod. A mod. B 


a 


46.3± 9.2 




41.8± 2.9 


- 


42.2 ±2.7 


- 


NR 


7.8 ± 7.8 




<3.5 


17 


< 3.5 




/o(980) 


6.2± 1.4 




4.1± 0.9 




4.8 ± 1.0 




/o(1370) 


2.3 ±1.7 




2.6 ±1.9 


3 


2.4 ± 1.3 




/o(1500) 






3.4 ±1.3 




3.4 ± 1.3 




(7r+7r~)s: 


(63 ± 12)* 


56.0 ± 3.9 (51.9 ±3.8)* 


108 


56 ±4 


50.4 ±2.7 56.6 ±5.2 


p(770) 


33.6±3.9 


30.8± 3.9 


20.0±2.5 


36 


25 ±4 


40 ±4 34.4 ±0.4 


/2(1270) 


19.4±2.5 


11.7± 1.9 


18.2±2.7 


21 


15.4 ± 2.5 


8.6±2.1 11.8 ±0.8 


p(1450) 


0.7±0.8 




<2.4 


1 


< 2.4 


1.5±0.6 0.8 ±0.4 




iKi.;-! 


OcS.T) 


90.1 


180 




100.2 lO.'-S.G 



[26] since they come from many different analyses. Our results are in general agreement with the 
others. The S'-wave is within the experimental results and the sum of the fit fractions of our models 
is very reasonable. 

The CLEO collaboration has performed an analysis where two models based on the mr scattering 
T-matrix were applied [14]. The model of Schechter (non isobar model, and chiral lagrangian with 
7 parameters) has similar fit fractions to our results although the coupling to the KK channel is 
not included in the model. The model of Achasov (isobar model with 12 parameters) produces 
results closest to [18] with large interference effects. The (7r''"7r'~)5 has a fit fraction of about 70% 
an the total sum is roughly 140%. 

Since we have performed an analysis that starts from the weak vertex, we can calculate the 
total branching ratio for the decay D'^ — tt+tt^tt"'" using the fit results for the parameters of the 
models. This gives 

model A : B{D+ 7r+7r-7r+) = (2.2^°;^) x 10-^, 
model B : B{D+ 7r+7r-7r+) = (2.51°;^) x 10-^ 

whose central values are smaller than the experimental average [26] 

^PDG(2P+ ^ 7r+7r-7r+) = (3.21 ± 0.19) x IQ-^^. (49) 
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The results are, however, coherent with the fact that we miss a part of the contribution in the a 
region. Concerning this discrepancy, one can advance some hypothesis. As aheady stated, this is 
Ukely due to the off-sheU effects that are not considered in our form factor. Another possible cause 
for the less prominent a peak and, consequently, for the smaller branching ratio, are three-body 
final state interactions. There are indications of three-body effects in K^ir^ir^ coming from 

a new analysis technique where the S'-wave is treated bin by bin in an almost model independent 
way [45, 46]. The Ktt S'-wave phase thus obtained is considerably different from the Kn scattering 
results and this has been interpreted as an indication of genuine three-body effects [47]. To our 
knowledge, this type of analysis has never been done for Z)+ — ^ tt'^tt^tt'^ but, in principle, three- 
body effects could also be important in this case. 

In conclusion, given the relative simplicity of our model, the small number of parameters and 
the fixed phases between the S'-wave and the the model is able to provide a fair description 
of the experimental data for Tr^vr^vr"^ within a unitary approach for the tttt final state 

interactions. Among many possible improvements to this work, one can think of introducing P 
wave form factors to describe the p(770)° and the p(1450)'^. 
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